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Abstract 

m 2 c 3 

The existence of electric fields close to their critical value E c — has been proved for massive 
cores of 10 7 up to 10 57 nucleons using a proton distribution of constant density and a sharp step 
function at its boundary [T] [2j [3] . We explore the modifications of this effect by considering a 
smoother density profile with a proton distribution fulfilling a Woods-Saxon dependence. The 
occurrence of a critical field has been confirmed. We discuss how the location of the maximum of 
the electric field as well as its magnitude is modified by the smoother distribution. 

1 Introduction 

One of the most active field of research has been to analyze a general approach to Neutron Stars 
based on the Thomas-Fermi ultrarelativistic equations amply adopted in the study of superheavy 
nuclei. The aim is to have a unified approach both to superheavy nuclei, up to atomic numbers 
of the order of 10 5 -10 6 , and to what we have called "Massive Nuclear Cores" . These cores are 
characterized by atomic number of the order of fO 57 , composed by neutrons, protons and electrons 
in /3-equilibrium, and expected to be kept at nuclear density by self gravity. 

The analysis of superheavy nuclei has historically represented a major field of research, guided 
by Prof. V. Popov and Prof. W. Greiner and their schools. This same problem was studied in 
the context of the relativistic Thomas-Fermi equation also by R. Ruffini and L.Stella, already in 
the 80s. A more recent approach [21 [3] has shown the possibility to extrapolate this treatment of 
superheavy nuclei to the case of Massive Nuclear Cores. 

The very unexpected result has been that also around these massive cores there is the distinct 
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possibility of having an electric field close to the critical value E c = e " ft , although localized in a 
very narrow shell of the order of the electron Compton wavelength. 

In all the mentioned works has been assumed a sharp profile for the proton distribution as 
given by a step function centered on the surface of the core; so modeling a sharp transition surface 
between for example, the core and the crust of neutron stars. In this work we model the transition 
surface in a smoother way by relaxing the sharp profile of the proton distribution and analyze the 
changes that it produce on the general properties of the system. 



2 The Relativistic Thomas— Fermi Equation 

Let us to introduce the proton distribution function f p (x) by mean of n p (x) = n P f p (x), where n p 
is the central number density of protons. We use the dimensionless unit x — (r — R c )/a, with 
a -1 = A-KOi\ e n pl A e is the electron Compton wavelength, R c the point where initial conditions 
are given (x = 0) and a is the fine structure constant. 

Using the Poisson's equation and the equilibrium condition for the gas of electrons 

Ep = m e c 2 ^/l + of — m e c 2 — eV = , (1) 
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Figure 1: Proton distribution function for 7 = 1.5, (3 « 0.0585749. 
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0.9662053 
0.97829293547 


-0.8680512263367902 
-0.899201 


1.38 x 10 36 
2.76 x 10 36 



Table 1 : Sets of initial conditions 



where e is the fundamental charge, x e the normalized electron Fermi momentum and V the 
electrostatic potential, we obtain the relativistic Thomas-Fermi equation 



Ce(x)- [ ^ X) 1]3/2 +/,(S)=0, 



C2) 



s l + R c /a J fi 

where n — 37r 2 A^rip and we have introduced the normalized electron chemical potential in 
absence of any field £ e = y/1 + x'j. For a given distribution function f p (x) and a central num- 
ber density of protons n^, the above equation can be integrated numerically with the boundary 
conditions 

£ e (0) = \/l + [/z5/ P (0)] 2/3 , ^(0)<0, (3) 

where 5 = n e (0)/n p (0). 

After integrating the TF equation, we can to calculate the neutron number density using the 
equilibrium condition of the direct and inverse (3 decay 



which results in 



e + p — > n + v , 



m n c^ n — m n c = m p c 2 £ p — m p c + eV 



(4) 



The electrostatic potential V is calculated using the equilibrium condition (TJJ . 

3 The Woods-Saxon— like Proton Distribution Function 

We propose a monotonically decreasing proton distribution function fulfilling a Woods-Saxon 
dependence 
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(5) 



7 + e^ x 

where 7 > and /3 > 0. In fig. [T]we show the proton distribution function for a particular set of 
parameters. 



4 Results of the Numerical Integration 

We have integrated numerically the eq.([2)| for several sets of parameters and initial conditions. 
To show the general behaviour of the system we show here the following two samples of initial 
conditions 

for the which we obtain respectively the physical quantities 
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Table 2: Physical quantities for the sets of parameters in table [T] 
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Figure 2: Electron and Proton Number Density for the sets of parameters in table [T] 

4.1 Number Densities 

In fig[2]we show the electron and proton number density together. We can see here that there is 
no local charge neutrality because the small difference between the two profiles. Nevertheless, this 
small difference creates a charge separation (see fig[3J) enough to produce huge electric fields (see 
figgj. 

We can also see how the system reaches indeed global charge neutrality in a very small scale 
as noted by Migdal et al. [I] in their classical paper. This scale have been calculated to be of the 
order of X n /y/a. 

4.2 The Charge Separation 

In order to see more clearly the difference between the electron and the proton profiles we have 
plotted in fig(3]the charge separation function given by 

A{x)= n p (x)-n e (x) (6) 

We see two well defined zones with opposite charge. In the first zone we have n p > n e so we 
have a positive charged shell while in the second one we have n v < n e and a negative charged shell 
develops. At the point n e = n p we have a maximum of the electric field, which is screened by the 
negative charged shell until reach global charge neutrality (see figUJ. 

4.3 The Electric Field 
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We have plotted in fig(4]the electric field in units of the critical field E c = — %— , namely 




Figure 3: Charge separation for the sets of parameters in table [T] 
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f = -T*«>- (7) 

We see that the electric field is overcritical but smaller in respect to the case of a sharp step 
proton distribution as used in [U|2][3]. Nevertheless, it is yet well above the critical field E c . 

The maximum of the electric field occur at the point where the transition from the positive 
charged shell to the negative charged one takes place, i.e., where n e — n v (see fig[3]). Of course 
because we have assumed a smoother proton profile we find also a smoother electric field about 
the maximum. We recall that in the case of sharp proton profile the first derivative of the electric 
field has a discontinuity on the point of charge inversion. 



5 Conclusions 

We confirm the existence of overcritical electric fields in the smooth transition surface of a massive 
nuclear core. The intensity of the electric field depends on the proton density mainly by two 
factors: the first one is the value of n p about the surface and the second one is how it changes 
about the surface (sharpness). In this line we note that the first factor depends strongly on a 
precise value of the so called 'melting density' and the correct value of the charge to mass ratio 
(Z/A) as given by (5— equilibrium, while for the second factor could be very important the surface 
tension as given for instance by the strong interaction. 
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